IN THE CLAIMS 

Please amend the claims as indicated: 1 . Canceled 



1 2. 


(currently amended) In a method for determining a fracture pressure gradient of a 


2 


subsurface region of earth formations comprising: 


3 


(a) 


obtaining seismic survey information about the subsurface region; 


4 


(b) 


identifying a plurality of interpreted seismic horizons of interest from the 


5 




obtained survey information; 


6 


(c) 


obtaining estimated seismic velocities corresponding to at least one 


7 




interval between at least one pair of said plurality of seismic horizons; 


8 


(d) 


calibrating the estimated seismic velocities to the parameter of interest 


9 


(e) 


using the results of said calibration and the obtained seismic velocities to 


10 




obtain said fracture pressure gradient at any location within the seismic 


11 




survey; 


12 


an improvement comprising displaying the parameter of interest on one of: 


13 


(0 


P- or S-wave velocity displays; 


14 


(ii) 


P-wave impedance displays; 


15 


(iii) 


S-wave impedance displays; 


16 


(iv) 


P-wave frequency attribute displays; 


17 


(v) 


S-wave frequency attribute displays; 


18 


(vi) 


P-wave phase attribute displays; 


19 


(vii) 


S-wave phase attribute displays; 


20 


(viii) 


density displays; 



2 



21 (ix) P-wave amplitude attribute displays; 

22 (x) S-wave amplitude attribute displays. 

3. Canceled 

4. Canceled. 

1 5 . (new) A computer system implementing a computer program comprising 

2 instructions for: 

3 (a) accessing subsurface seismic data that is at least one of: 

4 (A) a 2-D seismic line, and 

5 (B) a 3-D seismic volume, 

6 corresponding to a subsurface region; 

7 (b) displaying, editing and datuming a well log associated with the subsurface 

8 region and fitting a calibration curve to the log; 

9 (c) predicting fluid and rock pressures in the subsurface region based at least 

10 in part on the subsurface seismic data and results of the fitting; and 

1 1 (d) displaying results of the prediction. 
12 

1 6. (new) The computer system of claim 5 further comprising instructions for 

2 (i) identifying a plurality of seismic horizons from the subsurface seismic 

3 data; 

4 (ii) obtaining estimated seismic velocities corresponding to at least one 

5 interval between at least one pair of said plurality of seismic horizons; 

6 (iii) calibrating the estimated seismic velocities to the parameter of interest 



3 



7 

$ 7. (new) The computer system of claim 6 further comprising instructions for: 

2 calibrating the estimated seismic velocities using at least one of: 

3 (A) a function determined independently from the seismic data using regional 

4 information; 

5 (B) data from a well corresponding to a specific calibration location 

6 that is outside the areal extent of the seismic velocity data; and 

7 (C) data from a well corresponding to a specific calibration location within the 

8 areal extent of the seismic velocity data; and 

9 (D) data from a well corresponding to a specific calibration location within the 

10 areal extent of the seismic velocity data combined with the velocity data 

1 1 from the seismic survey for the same said location. 
12 

1 8. (new) The computer system of claim 7 further comprising instructions for 

2 calibrating the estimated seismic velocities based on estimation of an overburden- 

3 depth relationship that is determined by integrating density data from at least one 

4 of: 

5 I. a 1-D density function derived from density logs from at least one well; 

6 II. a 1-D density function derived from density data from drop cores, side- 

7 wall cores or conventional cores; 

8 m. a spatially varying density function based on the inversion of potential 

9 fields data, said potential fields data comprising at least one of gravity 
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10 data, and magnetic data; 

11 IV. a density function derived by performing inversion of at least one of PP 

12 data, and PS data. 
13 

14 9. (new) The computer system of claim 6 further comprising instructions for 

2 estimating the seismic velocities using at least one of: 

3 (A) Dip Move out (DMO), 

4 (B) Pre-stack time migration, 

5 (C) Pre-stack depth migration, 

6 (D) Multiple attenuation or suppression, 

7 (E) Reflection statics, 

8 (F) Refraction statics, 

9 (G) Wavefield reconstruction, and 

10 (H) combining undersampled seismic data from several gathers to form super- 

1 1 gathers. 
12 

1 1 0. (new) The computer system of claim 5 wherein the parameter of interest is 

2 selected from the group consisting of: (i) fluid pressure, (ii) fracture pressure, (iii) 

3 overburden pressure, (iv) effective stress, and, (v) excess pressure, defined as the 

4 difference between the actual fluid pressure and the normal pressure for the same 

5 depth. 
6 



5 



1 11. (new) The computer system of claim 1 0 further comprising instructions for 

2 performing further analysis at one of the following scales: 

3 A. the scale of a specific well location, defined as a projected one- 

4 dimensional path that will be used to place a wellbore into the subsurface, 

5 B. the scale of an exploration prospect, and 

6 C. the scale of a regional pressure evaluation to understand hydrocarbon 

7 systems. 
8 

9 12. (new) The computer system of claim 10 further comprising instructions for 

2 performing specific analyses of the subsurface including at least one of the 

3 following: 

4 A. basin modeling, 

5 B. hydrocarbon maturation in source rock intervals, 

6 C. lateral and vertical seal rock integrity, 

7 D. fault seal integrity, 

8 E. evaluation of fluid migration pathways in the subsurface, 

9 F. reservoir-specific lateral pressure changes, 

10 G. shallow water flow risk evaluation, and 

11 H. wellbore stability and failure analysis. 
12 

1 13. (new) The computer system of claim 5 wherein the seismic data is selected from 

2 the group consisting of: (i) P-wave land seismic data, (ii) S-wave land seismic 



6 



3 data, (iii) mode-converted S-wave land data (land PS data), (iv) mode-converted 

4 S-wave marine ocean-bottom cable data (marine PS data), (v) P-wave marine 

5 streamer data, and, (vi) P-wave marine ocean bottom cable data (PP data). 
6 

1 14. (new) The computer system of claim 6 further comprising instructions for 

2 estimating the seismic velocities from the group consisting of: 

3 (A) P-wave velocity data generated from normal moveout (NMO) velocity 

4 analysis, 

5 (B) P-wave or S-wave velocity data generated from coherency inversion 

6 analysis, 

7 (C) P-wave velocity generated from pre-stack inversion, 

8 (D) P-wave velocity generated from post-stack inversion, 

9 (E) S-wave velocity generated from pre-stack inversion, 

1 0 (F) S-wave velocity generated from post-stack inversion, 

1 1 (G) S-wave velocity data generated from normal moveout (NMO) velocity 

12 analysis, 

1 3 (H) P-wave velocity generated from tomography, 

14 (I) S-wave velocity generated from tomography, 

1 5 (J) P-wave velocity data from vertical seismic profiling (VSP), 

1 6 (K) P-wave velocity data from VSP look-ahead inversion, 

1 7 (L) S-wave velocity data from vertical seismic profiling (VSP), and 

1 8 (M) S-wave velocity data from VSP look-ahead inversion. 



7 



19 

20 15. (new) The computer system of claim 7 wherein the data from a well comprises at 

2 least one of: 

3 I. sonic logs, 

4 II. shear sonic logs, 

5 HI. density logs, 

6 IV. Lithology logs in the form of gamma ray or Spontaneous Potential logs, 

7 V. formation fluid pressure data, 

8 and 

9 VH. Leak off test data 
10 

1 16. (new) The computer system of claim 5 further comprising instructions for 

2 implementing at least one of: 

3 (i) the Eaton method, 

4 (ii) the effective stress method, and 

5 (iii) the equivalent depth method. 
6 

1 1 7. (new) The computer system of claim 6 further comprising instructions for 

2 calibrating the estimated seismic velocities using at least one of: 

3 (A) a linear curve fitting, 

4 (B) a curve fitting based on a power law, 

5 (C) a curve fitting based on exponentials, and 
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6 (D) a quadratic curve fitting. 
7 

S 18. (new) The computer system of claim 17 further comprising instructions for 

2 calibrating the estimated seismic velocities that enable use of an interactive 

3 display, the interactive display allowing modification of a displayed calibration 

4 curve and observing changes in a plurality of other displayed curves 

5 simultaneously, the interactive display further allowing at least one of: 

6 I. manually changing a coefficient or exponent in the fitting equation, and 

7 II. using a cursor to graphically manipulate the fitting curve. 
8 

1 19. (new) The computer system of claim 5 further comprising instructions for 

2 obtaining a porosity calibration curve using at least one of: 

3 (i) editing of erroneous porosity values for at least one well; 

4 (ii) displaying in a single plot, a depth-correlated porosity and velocity data 

5 for the at least one well; 

6 (iii) curve fitting of the velocity and porosity data using at least one of: 

7 La linear curve fitting, 

8 n. a curve fitting based on a power law, 

9 HI. a curve fitting based on exponentials, and 
10 IV. a quadratic curve fitting. 

11 



1 20. (new) The computer system of claim 5 further comprising instructions fo displaying 



9 



2 the parameter of interest on one of: 

3 (i) stacked seismic section, and 

4 (ii) migrated seismic section. 
5 

6 21 . (new) The computer system of claim 5 further comprising instructions for enabling 

2 display of the parameter of interest interactively and simultaneously on at least 

3 one of: 

4 (i) a seismic display; 

5 (ii) a velocity versus depth display including a velocity function for a specific 

6 location and a calibration function for velocity versus effective stress, 

7 (iii) a stress versus depth display including the overburden stress calibration 

8 for said specific location and the effective stress calculated from the 

9 velocity versus depth display, 

10 (iv) a pressure-gradient versus depth display including the fracture gradient or 

1 1 overburden gradient, the fluid pressure gradient calculated from the stress 

12 versus depth display, and pressure data points from a well that applies to 

13 said specific location. 
14 

1$ 22. (new) The computer system of claim 5 further comprising instructions enabling at 

2 least one of: 

3 (i) deletion or correction for zones of abnormal velocity caused by the 

4 presence of hydrocarbon fluids; and 
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5 (ii) deletion of zones of abnormal velocity caused by the presence of 

6 nonclastic rocks. 
7 

$ 23. (new) The computer system of claim 5 further comprising instructions for at least one 

2 of: 

3 (i) determining an overburden stress from a density function based on 

4 inversion of 2-D or 3-D potential fields; 

5 (ii) determining an overburden stress from a density function derived by 

6 performing a simultaneous inversion of PP and PS data; 

7 (iii) correcting for zones of abnormal density caused by the presence of 

8 hydrocarbon fluids by inserting a correct density for said zones and 

9 recalculating the overburden at the specific location; 

1 0 (iv) correcting for zones of abnormal velocity caused by the presence of 

1 1 hydrocarbon fluids by deleting said zones from the calculation of the 

1 2 parameter of interest; 

1 3 (v) correcting for zones of abnormal density caused by the presence of non- 

14 clastic rocks by inserting a correct density for said zones and recalculating 

15 the overburden at the specific location; and 

16 (vi) correcting for zones of abnormal velocity caused by the presence of non- 

17 clastic rocks by deleting said zones from the calculation of the parameter 

18 of interest. 



19 
20 

11 



1 24. (new) The computer system of claim 6 further comprising instructions for estimation 

2 of an overburden-depth relationship that is determined by integrating density data 

3 obtained by inversion of 2-D or 3-D potential fields data. 
4 

5 25. (new) The computer system of claim 6 further comprising instructions for estimation 

2 of an overburden-depth relationship by integrating density data obtained by 

3 inversion of at least one of 2-D or 3-D seismic data, and wherein said seismic data 

4 further comprise at least one of PP data and PS data. 
5 

1 26. (new) A machine readable medium comprising instructions for: 

2 (a) accessing subsurface seismic data that is at least one of: 

3 (A) a 2-D seismic line, and 

4 (B) a 3-D seismic volume, 

5 corresponding to a subsurface region; 

6 (b) displaying, editing and datuming a well log associated with the subsurface 

7 region and fitting a calibration curve to the log; 

8 (c) predicting fluid and rock pressures in the subsurface region based at least 

9 in part on the subsurface seismic data and results of the fitting; and 
1 0 (d) displaying results of the prediction. 

11 

1 27. (new) The machine readable medium of claim 26 further comprising instructions 

2 for 
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3 (i) identifying a plurality of seismic horizons from the subsurface seismic 

4 data; 

5 (ii) obtaining estimated seismic velocities corresponding to at least one 

6 interval between at least one pair of said plurality of seismic horizons; 

7 (iii) calibrating the estimated seismic velocities to the parameter of interest 
8 

9 28. (new) The machine readable medium of claim 27 further comprising instructions 

2 for: 

3 calibrating the estimated seismic velocities using at least one of: 

4 (A) a function determined independently from the seismic data using regional 

5 information; 

6 (B) data from a well corresponding to a specific calibration location 

7 that is outside the areal extent of the seismic velocity data; and 

8 (C) data from a well corresponding to a specific calibration location within the 

9 areal extent of the seismic velocity data; and 

10 (D) data from a well corresponding to a specific calibration location within the 

1 1 areal extent of the seismic velocity data combined with the velocity data 

12 from the seismic survey for the same said location. 
13 

1 29. (new) The machine readable medium of claim 28 further comprising instructions 

2 for calibrating the estimated seismic velocities based on estimation of an 

3 overburden-depth relationship that is determined by integrating density data from 
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4 


at least one of: 


5 


I. 


a 1-D density function derived from density logs from at least one well; 


6 


n. 


a 1-D density function derived from density data from drop cores, side- 


7 




wall cores or conventional cores; 


8 


m. 


a spatially varying density function based on the inversion of potential 


9 




fields data, said potential fields data comprising at least one of gravity 


10 




data, and magnetic data; 


11 


rv. 


a density function derived by performing inversion of at least one of PP 


12 




data, and PS data. 


13 






14 30. 


(new) The machine readable medium of claim 27 further comprising instructions 


2 


for estimating the seismic velocities using at least one of: 


3 


(A) 


Dip Move out (DMO), 


4 


(B) 


Pre-stack time migration, 


5 


(C) 


Pre-stack depth migration, 


6 


(D) 


Multiple attenuation or suppression, 


7 


(E) 


Reflection statics, 


8 


(F) 


Refraction statics, 


9 


(G) 


Wavefield reconstruction, and 


10 


(H) 


combining undersampled seismic data from several gathers to form super- 


11 




gathers. 


12 







1 31. (new) The machine readable medium of claim 26 wherein the parameter of 

2 interest is selected from the group consisting of: (i) fluid pressure, (ii) fracture 

3 gradient, (iii) overburden pressure, (iv) effective stress, and, (v) excess pressure, 

4 defined as the difference between the actual fluid pressure and the normal 

5 pressure for the same depth. 
6 

1 32. (new) The machine readable medium of claim 3 1 further comprising instructions 

2 for performing further analysis at one of the following scales: 

3 A. the scale of a specific well location, defined as a projected one- 

4 dimensional path that will be used to place a wellbore into the subsurface, 

5 B. the scale of an exploration prospect, and 

6 C. the scale of a regional pressure evaluation to understand hydrocarbon 

7 systems. 
8 

$ 33. (new) The machine readable medium of claim 31 further comprising instructions for 

2 performing specific analyses of the subsurface including at least one of the 

3 following: 

4 A. basin modeling, 

5 B. hydrocarbon maturation in source rock intervals, 

6 C. lateral and vertical seal rock integrity, 

7 D. fault seal integrity, 

8 E. evaluation of fluid migration pathways in the subsurface, 
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9 F. reservoir-specific lateral pressure changes, 

10 G. shallow water flow risk evaluation, and 

1 1 EL wellbore stability and failure analysis. 
12 

1 34. (new) The machine readable medium of claim 26 wherein the seismic data is 

2 selected from the group consisting of: (i) P-wave land seismic data, (ii) S-wave 

3 land seismic data, (iii) mode-converted S-wave land data (land PS data), (iv) 

4 mode-converted S-wave marine ocean-bottom cable data (marine PS data), (v) P- 

5 wave marine streamer data, and, (vi) P-wave marine ocean bottom cable data (PP 

6 data). 
7 

1 35. (new) The machine readable medium of claim 27 further comprising instructions 

2 estimating the seismic velocities from the group consisting of: 

3 (A) P-wave velocity data generated from normal moveout (NMO) velocity 

4 analysis, 

5 (B) P-wave or S-wave velocity data generated from coherency inversion 

6 analysis, 

7 (C) P-wave velocity generated from pre-stack inversion, 

8 (D) P-wave velocity generated from post-stack inversion, 

9 (E) S-wave velocity generated from pre-stack inversion, 

1 0 (F) S-wave velocity generated from post-stack inversion, 

1 1 (G) S-wave velocity data generated from normal moveout (NMO) velocity 
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analysis, 

(H) P-wave velocity generated from tomography, 

(I) S-wave velocity generated from tomography, 

(J) P-wave velocity data from vertical seismic profiling (VSP), 

(K) P-wave velocity data from VSP look-ahead inversion, 

(L) S-wave velocity data from vertical seismic profiling (VSP), and 

(M) S-wave velocity data from VSP look-ahead inversion. 

(new) The machine readable medium of claim 28 wherein the data from a well 
comprises at least one of: 

I. sonic logs, 

II. shear sonic logs, 
HI. density logs, 

IV. Lithology logs in the form of gamma ray or Spontaneous Potential logs, 
and 

V. formation fluid pressure data, 

(new) The machine readable medium of claim 26 further comprising instructions 
for implementing at least one of: 

(i) the Eaton method, 

(ii) the effective stress method, and 

(iii) the equivalent depth method. 

17 



6 

1 38. (new) The machine readable medium of claim 27 further comprising instructions 

2 for calibrating the estimated seismic velocities using at least one of: 

3 (A) a linear curve fitting, 

4 (B) a curve fitting based on a power law, 

5 (C) a curve fitting based on exponentials, and 

6 (D) a quadratic curve fitting. 
7 

$ 39. (new) The machine readable medium of claim 38 further comprising instructions 

2 for calibrating the estimated seismic velocities that enable use of an interactive 

3 display, the interactive display allowing modification of a displayed calibration 

4 curve and observing changes in a plurality of other displayed curves 

5 simultaneously, the interactive display further allowing at least one of: 

6 I. manually changing a coefficient or exponent in the fitting equation, and 

7 II. using a cursor to graphically manipulate the fitting curve. 
8 

$ 40. (new) machine readable medium of claim 26 further comprising instructions for 

2 obtaining a porosity calibration curve using at least one of: 

3 (i) editing of erroneous porosity values for at least one well; 

4 (ii) displaying in a single plot, a depth-correlated porosity and velocity data 

5 for the at least one well; 

6 (iii) curve fitting of the velocity and porosity data using at least one of: 
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7 La linear curve fitting, 

8 n. a curve fitting based on a power law, 

9 HI. a curve fitting based on exponentials, and 
10 IV. a quadratic curve fitting. 

11 

12 41 . (new) The machine readable medium of claim 26 further comprising instructions fo 

2 displaying the parameter of interest on one of: 

3 (i) stacked seismic section, and 

4 (ii) migrated seismic section. 
5 

6 42. (new) The machine readable medium of claim 26 further comprising instructions for 

2 enabling display of the parameter of interest interactively and simultaneously on 

3 at least one of: 

4 (i) a seismic display; 

5 (ii) a velocity versus depth display including a velocity function for a specific 

6 location and a calibration function for velocity versus effective stress, 

7 (iii) a stress versus depth display including the overburden stress calibration 

8 for said specific location and the effective stress calculated from the 

9 velocity versus depth display, 

1 0 (iv) a pressure-gradient versus depth display including the fracture gradient or 

1 1 overburden gradient, the fluid pressure gradient calculated from the stress 

12 versus depth display, and pressure data points from a well that applies to 
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13 said specific location. 

14 

IS 43. (new) The machine readable medium of claim 6 further comprising instructions 

2 enabling at least one of: 

3 (i) deletion or correction for zones of abnormal velocity caused by the 

4 presence of hydrocarbon fluids; and 

5 (ii) deletion of zones of abnormal velocity caused by the presence of 

6 nonclastic rocks. 
7 

$ 44. (new) The machine readable medium of claim 26 further comprising instructions for 

2 at least one of: 

3 (i) determining an overburden stress from a density function based on 

4 inversion of 2-D or 3-D potential fields; 

5 (ii) determining an overburden stress from a density function derived by 

6 performing a simultaneous inversion of PP and PS data; 

7 (iii) correcting for zones of abnormal density caused by the presence of 

8 hydrocarbon fluids by inserting a correct density for said zones and 

9 recalculating the overburden at the specific location; 

1 0 (iv) correcting for zones of abnormal velocity caused by the presence of 

1 1 hydrocarbon fluids by deleting said zones from the calculation of the 

1 2 parameter of interest; 

13 (v) correcting for zones of abnormal density caused by the presence of non- 
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14 clastic rocks by inserting a correct density for said zones and recalculating 

15 the overburden at the specific location; and 

1 6 (vi) correcting for zones of abnormal velocity caused by the presence of non- 

17 clastic rocks by deleting said zones from the calculation of the parameter 

18 of interest. 
19 

20 45. (new) The machine readable medium of claim 27 further comprising instructions for 

2 estimation of an overburden-depth relationship that is determined by integrating 

3 density data obtained by inversion of 2-D or 3-D potential fields data. 
4 

$ 46. (new) The machine readable medium of claim 27 further comprising instructions for 

2 estimation of an overburden-depth relationship by integrating density data 

3 obtained by inversion of at least one of 2-D or 3-D seismic data, and wherein said 

4 seismic data further comprise at least one of PP data and PS data. 
5 

1 47. (new) A computer program executed by a computer for: 

2 (a) accessing subsurface seismic data that is at least one of: 

3 (A) a 2-D seismic line, and 

4 (B) a 3-D seismic volume, 

5 corresponding to a subsurface region; 

6 (b) displaying, editing and datuming a well log associated with the subsurface 

7 region and fitting a calibration curve to the log; 
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8 (c) predicting fluid and rock pressures in the subsurface region based at least 

9 in part on the subsurface seismic data and results of the fitting; and 

1 0 (d) displaying results of the prediction on a display device associated with the 

1 1 computer. 
12 

1 48. (new) The computer program of claim 47 further comprising instructions for: 

2 (i) identifying a plurality of seismic horizons from the subsurface seismic 

3 data; 

4 (ii) obtaining estimated seismic velocities corresponding to at least one 

5 interval between at least one pair of said plurality of seismic horizons; 

6 (iii) calibrating the estimated seismic velocities to the parameter of interest 
7 

$ 49. (new) The computer program of claim 48 further comprising instructions for: 

2 calibrating the estimated seismic velocities using at least one of: 

3 (A) a function determined independently from the seismic data using regional 

4 information; 

5 (B) data from a well corresponding to a specific calibration location 

6 that is outside the areal extent of the seismic velocity data; and 

7 (C) data from a well corresponding to a specific calibration location within the 

8 areal extent of the seismic velocity data; and 

9 (D) data from a well corresponding to a specific calibration location within the 

10 areal extent of the seismic velocity data combined with the velocity data 
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1 1 from the seismic survey for the same said location. 

12 

13 



23 



